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Progress to Date 


This report documents the progress that has been made in the proposed research work over the 
past year. The goals for year one have been completed. In addition, substantial progress has 

been made toward fulfilling the goals of year two. Work has also begun on additional areas that 
have been identified. 


The work that has been completed to date is: 


Completed before semi-annual report 

1. A simplified rotor model with a flexible shaft and backup bearings has been developed 
The model is based upon the work of Ishii and Kirk. Parameter studies of the behavior 
of this model are currently being conducted. 

2. A simple rotor model which includes a flexible disk and bearings with clearance has been 
developed and the dynamics of the model investigated. The study consists of simulation 
work coupled with experimental verification. The work is documented in the attached 
paper. A copy of the paper is included in the Appendix. 

T A rotor model based upon the T-501 engine has been developed which includes backup 
bearing effects. The dynamics of this model are currently being studied with the 
objective of verifying the conclusions obtained from the simpler models. 

1. Parallel simulation runs are being conducted using an ANSYS based finite element model 
ot the T-501 developed by Dr. Charles Lawrence. 


Progress since semi-annual report 

1. The magnetic bearing test rig is currently floating and dynamics/control tests are 
currently being conducted. Photographs of the various test rigs that have been assembled 
for experimental testing for this project are shown in Figures 1 - 3. 

2. A paper has been written that documents the work using the T-501 engine model. The 
paper has been accepted for presentation at the Symposium on Nonlinear and Stochastic 
Dynamics, to be held at the 1994 ASME Winter Annual Meeting, November 13-18, 

1994, Chicago, Illinois. A copy of the paper is included in the Appendix. 

3. Work has continued with the simplified model. A paper that fully documents this work 
is currently being written and will be submitted to the ASME Turbo Expo June 5-8 

1995, Houston, Texas. ’ 
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4. The finite element model is currently being modified to include the effects of foundation 
dynamics. 

5. A literature search for material on foil bearings has been conducted. 

6. A finite element model is being developed for a magnetic bearing in series with a foil 
backup bearing. 
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Figure 3: Transient Drop Rotordynamics Test Rig 
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Abstract 

This study presents an investigation of the 
dynamics of a rotor system with bearing clearance. 
Of particular interest is the influence of such effects 
on coupled shaft/disk vibration. Experimental 
results for a rotor system with a flexible disk 
are presented and compared to predictions from a 
simulation model. Some insights and conclusions 
are obtained with regard to the conditions under 
which such vibration may be significant. 

Nomenclature 

e = imbalance eccentricity 
fni, f n 2 = nonlinear bearing forces 
r i — gyroscopic mass influence ratio 
To = disk mass influence ratio 
Xi , xo ~ shaft degrees of freedom 
X 3 , x,[ = disk degrees of freedom 
A = bearing clearance 
f 2 = rotor speed 

= natural frequency of rotor support in Xi 
direction 

u>o = natural frequency of rotor support in X 2 
direction 

= natural frequency of rotor disk 

T — fit 

£1 = damping ratio of rotor support in Xi 
direction 

£2 = damping ratio of rotor support in X 2 
direction 

£3 = damping ratio of rotor disk 



Background and Motivation 

Rotor systems typically consist of a shaft with 
one or more bladed disks attached. Disk and blade 
vibration are issues of important concern for stress 


analysts. Complex finite element models are devel- 
oped to assess the dynamic stresses in such compo- 
nents and insure the integrity of the design. How- 
ever, models for rotordynamical analyses are typ- 
ically developed assuming that disk flexibility '’ef- 
fects are negligible. It is generally presumed that 
the disk is sufficiently rigid so as to not signifi- 
cantly impact rotor vibration over the operating 
region. 

There is a fairly large body of work documented 
in the literature concerning studies of disk flexi- 
bility on rotors and turbomachinery. A discussion 
of this previous work with regard to linear rotor 
systems is presented in Flowers and Ryan (1991). 
An excellent source for a comprehensive review of 
work in this area is presented by Davis (1989). 

Previous investigators into the area of coupled 
rotor/disk vibration have noted that disk flexi- 
bility has little effect on critical speeds but that 
it may significantly influence higher natural fre- 
quencies of the rotor system [Chivens and Nel- 
son (1975)]. One can draw the conclusion from 
these results that synchronous vibration due to 
imbalance will be little affected by disk flexibil- 
ity. However, there are sources of higher fre- 
quency excitation that could serve to excite these 
natural frequencies. Perhaps the most obvious 
are multi-synchronous effects corresponding to the 
blade pass frequency (from fluid forces impinging 
on the rotor blades). Another potential source is 
from nonlinearities that may be present in the sys- 
tem. For example, nonlinear bearing forces due to 
clearance effects may result in supersvnchronous 
rotor vibration. There are quite a number of stud- 
ies in the literature concerning the effects of bear- 
ing clearance (and the related phenomena of rub- 
bing) on rotordynamical behavior. Some of the 
works that have most influenced the current study 
were conducted by Johnson (1962), Black (1968), 
Ehrich (1966), and Childs (1978), and Muszvnska 
(1984). 



Advanced designs for many types of high speed 
rotating machinery that use magnetic bearings for 
support have been proposed and are in the de- 
velopment and construction stages. There are a 
number of such systems already in commercial 
use. Rotor systems supported by magnetic bear- 
ings must have backup bearings to provide support 
under overload conditions or if the magnetic bear- 
ing fails. Backup bearings are characterized by a 
clearance between the rotor and the bearing such 
that contact does not occur under normal operat- 
ing conditions. As a result, issues related to the 
effects of bearing clearance on rotordynamical be- 
havior are of current concern. The objective of 
the present work is to develop an understanding 
of possible coupling between the dynamics of disk 
and shaft that may be induced by bearing clear- 
ance effects and to provide guidance to designers 
concerned with such systems. 

Experimental Model 

In order to investigate whether bearing clear- 
ance can lead to significant coupling between rotor 
and disk vibration, experimental tests were per- 
formed with a rotor test rig. A drawing of the test 
rig is shown in Figure 1 . 



Horizontal Sensors 


Figure 1 Experimental Test Apparatus 

The rotor used in this study has two basic 
components: a flexible disk and a shaft. The steel 
shaft is 0.375 in. diameter and 18 inches long. The 
flexible disk is a circular aluminum plate 0.0125 
inches thick and 14.0 inches in diameter. The 
natural frequency of its lowest one nodal diameter 
bending mode is about 40 Hz. It is attached with 
bolts to a 2.0 inch diameter hub that is fixed on at 



Figure 2 Shaft Right End Support Device 


the center of a 0.375 inch diameter steel shaft 16.0 
inches in length. 

The rotor shaft is supported by bushings at two 
ends. The left end of the shaft is placed directly in 
the bushing base. This support provides a force 
but only a minimal couple so that the shaft is 
effectively free to pivot about this point. The right 
bearing is supported by a special device, which 
is designed to provide stiffness and to simulate 
bearing clearance. A diagram for this device is 
shown in figure 2. It has two sets of springs. 
The softer springs are to support the weight of 
the rotor, and the harder springs act as the main 
stiffness for the system. The stiffness in the x and 
y directions of each spring set may be different. 
The clearance is adjustable by turning the inner 
tapered ring in or out. In this experiment, the 
spring constant of the soft spring is 1.9 lb/in. The 
spring constant of the hard spring is 8.1 lb /in. in 
the xi (horizontal) direction and 19.7 lb/in. in the 
X 2 (vertical) direction. 

The rotor is driven by a speed adjustable motor. 
The speed of the rotor can be controlled by turning 
the speed control knob on the front panel of 
the control box. The measurement and analysis 
system comprises proximity displacement sensors, 
sensor conditioners and a signal analyzer. A 
vertical and a horizontal displacement sensor were 
mounted to pick up the displacement signals to 
form orbit plots. A third transducer was used 








to sense the disk vibration. During the test, 
displacement signals of the shaft were sent to the 
signal analyzer, where the orbit trajectories were 
recorded and the frequency components of the 
vibration were analyzed. 


Simulation Model 

The simulation model is similar to that devel- 
oped in an earlier study [Flowers and Wu (1993)]. 
The primary differences are that rotational stiffen- 
ing of the disk have been taken into account and 
that damping is included in the bearing support 
forces. A schematic diagram of the model is shown 
in figure 3. It consists of a rigid shaft, a rigid hub, 
a flexible disk and a support with a symmetric 
clearance. The equation development is based on 
the following considerations: 

(1) The disk is assumed to flex only in the lateral 
direction, 

(2) Only rotational vibration is considered. This is 
because only the rotational motion is coupled 
with the lateral disk vibration. One nodal 
diameter disk vibration is assumed. 

(3) The rotor speed is constant. 


Flexible Disk Clearance Mechanism 



Figure 3 Simulation Model 

After some mathematical manipulations, the 
dimensionless equations of motion for the whole 
system can be obtained as 

*i + r x x 2 + r 2 x 3 + 2 x 4 + / nl = cos r (1) 

x 2 - tix x + r 2 x 4 - 2 x 3 + f n2 ~ sin r (2) 


x 3 + 2 x 4 + 2^3^32:3 + (u>3 — 1)2:3 

= ~(^i + 2 ^ 2 ) (3) 


x 4 — 2 x 3 + 2 £ 3 uj 3 x 4 + (u>3 - l)x 4 

= -(x 2 - 2 i 1 ) ( 4 ) 


fnl 



In 2 = —<t> 



e 



“7=t=7=t)( w i x i ) + 26«ixi 

V x i + x 2 J 

7 = T , 0 (^2 * 2 ) + 2<f 2 W 2 X 2 

V X 1 + x 2 


4 >= 1 + xj > 6 

0 otherwise 

In the above equations, x lt x 2 , x 3 and x 4 are 
scaled by e, the imbalance eccentricity. Time is 
scaled by the rotor speed, S7. f nl and / n2 are the 
nonlinear moments due to contacting between the 
rotor and the bearings. 

The effective stiffness of a spinning disk can be 
strongly influenced by the spin speed. Based upon 
earlier work [Wu and Flowers (1992)], the lowest 
natural frequency of such a disk can be written as 
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19 + i/ 
16 
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where u 3 q is the natural frequency of the nonspin- 
ning disk. This expression is used in equations 3 
and 4 to account for spin stiffening effects. 


Comparison of Experiment and Simulation 

A number of tests were conducted using the ex- 
perimental setup described in the previous section. 
In addition, studies for differing parametric con- 
figurations have been conducted using the simu- 
lation model described above. The basic results 
are very similar to those obtained by Flowers and 
Wu (1993). The results presented here are typical. 
The following discussion is directed at comparing 
the predictions of the simulation model with exper- 
imentally observed responses, with the objective of 
obtaining insight into the behavior of flexible disk 
rotor systems with bearing clearance effects. 

First, the linear characteristics of the system 
were examined. Figure 4 shows experimentally 
obtained data of the natural frequencies for the 
coupled disk/shaft vibration mode as a function 
of rotor speed. Using the data from such tests 
and additional measurements and calculations, the 
linear stiffness, damping, and mass characteristics 
of the experimental rotor system were identified. 
The numerical values are shown in Table 1. 




Parameter 

Value 

r i 

0.5 

*2 

0.1 

6 

0.06 

£2 

0.06 

& 

0.01 

U>i 

0.55 

U> 2 

0.85 1 

U>3 

1.8 

V 

0.3 

s 

1.2 


Table 1 Simulation Model Parameters 


Figure 4 also shows the natural frequencies 
as a function of rotor speed obtained from the 
simulation model. These results agree relatively 
well with the experimentally obtained data and 
serve to validate the structural parameters selected 
for the simulation model. ft, 3ft, and 5ft lines are 
also shown on the figure. The intersections of these 
lines with the natural frequency curve provides 
significant insight into the nonlinear behavior of 
the flexible disk rotor system, as is discussed 
below. 



Figure 4 Rotor System Natural Frequencies 


Next, the nonlinear behavior of the system was 
investigated. A clearance value was selected and 
the rotor speed was slowly increased over a range 
of rotor speeds, with the resulting response ampli- 
tudes and frequencies measured and recorded. 

It is important to note that many vibration fre- 
quencies are possible. During the course of the 
study, lft - 6 ft frequency components were ob- 
served for certain parametric configurations. In 
addition, frequency components that are fractions 
of the rotor speed were observed for certain speed 
ranges and initial conditions. However, for a ro- 
tor with a symmetric clearance and minimal fric- 
tional effects, the odd integer multiples of the ro- 
tor speed appear to generally be the primary su- 
persynchronous frequency components. The cur- 
rent discussion will concentrate on those frequen- 
cies. However, the basic conclusions should be di- 
rectly applicable to other supersynchronous vibra- 
tion frequencies as well. 

Figure 5 shows the experimentally measured 
amplitudes of the ft, 3ft, and 5ft components of the 
shaft response for various rotor speeds. Simulation 
results are also shown on figure 5. These results 
were obtained using the harmonic balance method 
and were verified at selected points through di- 
rect numerical integration of the equations of mo- 
tion. There is relatively good agreement between 
the simulation results and the corresponding ex- 
perimental data. Certainly the basic trends match 
quite well. 

An important observation that can be made 
from examination of these figures is that the peaks 
of the respective components correspond to the 
the intersections of the ft, 3ft, and 5ft lines with 
the natural frequency curves. It appears that the 
nonlinear effects are serving to excite the coupled 
modes of the rotor/disk system. This is true for 
both the forward and backward whirl modes of the 
rotor system. 

The peaks occurring at a rotor speed of 0.51 cor- 
respond to the intersections of the 3ft line with the 
natural frequency curve. They relate to a back- 
ward whirl mode. The peaks occurring at rotor 
speeds of 0.33 and 0.485 correspond to the intersec- 
tions of the 5ft line with the the natural frequency 
curve. Those occurring at 0.33 relate to a back- 
ward whirl mode and those occurring at 0.485 re- 
late to a forward whirl mode. Note that the ampli- 
tudes of the supersynchronous components associ- 
ated with the forward whirl mode of the disk/shaft 
vibration are much lower than those associated 
with the backward whirl mode. This result was 
observed for all the parametric configurations ex- 
amined in the course of this study. 

The ft, 3ft, and 5ft components contribute 
(more or less significantly) to the overall response 
at the corresponding intersection point. However, 
it is important to note that the ft component is 




increased rather dramatically at the intersection 
points as a result of the presence of disk flexibility. 
Apparently, the nonlinear coupling between the 
various frequency components serves to produce 
this effect. 


Simulation Model 



Figure 5. a Component of Shaft Vibration 
Response 


Conclusions 

A study of disk/shaft vibration induced by bear- 
ing clearance effects has been presented. Both ex- 
perimental tests and simulation results have been 
presented. The responses predicted by the sim- 
ulation model and those observed experimentally 
agree quite well. The behavior of both the ex- 
perimental test rig and the simulation model was 
quite sensitive to changes in parametric configu- 
ration. As a result it is quite difficult to predict 
exactly how a certain system is going to behave 
in the presence of a bearing clearance effect, i.e., 
whether or not coupled disk/shaft vibration will 
be significantly excited. However, a few general 
conclusions and guidelines can be drawn from this 
study. 

L An understanding of the mechanism for cou- 
pled disk/shaft vibration induced by clearance 
effects has been obtained. The nonlinear ef- 
fects have served to produce superharmonics 
that excite coupled disk shaft modes. Con- 
versely, the additional degree of freedom pro- 


Simulation Model 



Figure 5*b Component of Shaft Vibration 
Response 


Simulation Model 



Figure 5.c 5fi Component of Shaft Vibration 
Response 



vided by the disk flexibility has served to ex- 
aggerate the frequency components of the re- 
sponse, resulting in higher amplitudes at the ro- 
tor speeds corresponding to the respective nat- 
ural frequency /multi-synchronous line intersec- 
tion points. 

2. Bearing clearance effects primarily serve to ex- 
cite backward whirl modes for coupled disk/shaft 
vibration. Forward whirl modes may be ex- 
cited but the associated amplitudes appear to 
be much lower. 

3. The bearing support stiffnesses used in this 
study were not symmetric. These nonsymme- 
tries, together with the clearance, have served 
to produce the supersynchronous vibration that 
excites the coupled disk/shaft vibration modes. 

4. It appears to be relatively difficult to excite cou- 
pled disk/shaft vibration with bearing clearance 
effects. The occurrence of such behavior was 
very sensitive to rotor speed, as judged by the 
relative difficulty encountered during the exper- 
imental work in obtaining the peak amplitude 
responses. This conclusion is verified by the sim- 
ulation work that indicates that the behavior 
occurs only over a very limited ranges of rotor 
speed. The supersynchronous frequency must 
almost exactly coincide with a natural frequency 
in order to excite the behavior. 

5. From a design perspective, the Campbell dia- 
gram is a useful tool to predict when such be- 
havior may occur. The intersections of the cor- 
responding supersynchronous line with the cou- 
pled disk/shaft vibration mode frequency curve 
will indicate at what rotor speeds peak responses 
for that frequency component are likely to occur. 
Whether or not such behavior will actually oc- 
cur and the relative significance of such effects 
depends very strongly on the imbalance, clear- 
ance, damping, and stiffness values. 
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ABSTRACT 

This paper investigates the steady-state responses of a ro- 
tor system supported by auxiliary bearings in which there 
is a clearance between the rotor and the inner race of the 
bearing. A simulation model based upon the rotor of a pro- 
duction jet engine is developed and its steady-state behavior 
is explored over a wide range of operating conditions for var- 
ious parametric configurations. Specifically, the influence of 
rotor imbalance, support stiffness and damping is studied. 
It is found that imbalance may change the rotor responses 
dramatically in terms of frequency contents at certain oper- 
ating speeds. Subharmonic responses of 2nd order through 
10th order are all observed except the 9th order. Chaotic 
phenomenon is also observed. Jump phenomena (or double- 
valued responses) of both hard— spring type and soft— spring 
type are shown to occur at low operating speeds for sys- 
tems with low auxiliary bearing damping or large clearance 
even with relatively small imbalance. The effect of friction 
between the shaft and the inner race of the bearing is also 
discussed. 

NOMENCLATURE 

Cb = auxiliary bearing support damping, lb.s 2 /in. 

Cb^p = auxiliary bearing torsional damping, lb.in.s 
F n = normal force, lb 
F t = friction force, lb 

F* = external force vector acting on the rotor in X direction 
Fy — external force vector acting on the rotor in Y direction 
I«a = rotor inertia matrix 


Jb = moment of inertia of auxiliary bearing, Ib.in^ 2 
Kb = auxiliary bearing support stiffness, lb/in. 

Kc = contact stiffness, lb/in. 

Mb = auxiliary bearing mass, lb.s 2 /in. 

Mk = mass of kth rotor element, lb.s 2 /in. 

N = total number of modes considered 
MB 1 = node number at auxiliary bearing 
N B 2 = node number at auxiliary bearing #2 
Qx = rotor modal coordinate vector in X direction 
Qy = rotor modal coordinate vector in Y direction 
Rb = radius of auxiliary bearing bore, in. 

Rm = radius of auxiliary bearing pitch, in. 

Rr = radius of rotor journal, in. 

X/i = rotor physical coordinate vector in X direction 
Y/j = rotor physical coordinate vector in Y direction 
e = rotor imbalance eccentricity, in. 
g = gravitational acceleration, in./s 2 
t = time, s 

A = deformation at the contact point, in. 

r = * r i a * 

® = rotor free-free modal rotation matrix 
Q = rotor operating speed, rad/s 
$ = rotor free-free modal displacement matrix 
S = Rb — Rr, auxiliary bearing clearance, in. 

/i = dynamic friction coefficient 
= rolling friction coefficient 

xps — angular displacement of auxiliary bearing inner-race 
C = modal damping coefficient 
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INTRODUCTION 

One of the most innovative developments in the turbo- 
machinery field involves the use of active magnetic bearings 
(AMB) for rotor support. This technology provides the po- 
tential for significant improvements in the dynamic behavior 
of rotor systems, allowing for loading, eccentricity, shaft po- 
sition and vibration to be continuously monitored and con- 
trolled. In order to protect the soft iron cores of the magnetic 
bearings and to provide rotor support in the event of failure 
of the bearing or during an overload situation, backup (or 
auxiliary) bearings, with a clearance between the rotor and 
the inner race of the bearing, are usually included in the ro- 
tor design. This clearance introduces a nonlinear dynamical 
feature which may significantly impact the behavior of the 
rotor. 

Magnetic bearing systems appear to provide particularly 
great promise for use in aerospace applications. There are 
active programs at many of the major jet engine manufac- 
turers to develop engines supported by magnetic bearings. 
Safety is a major concerns in any aeronautical design. To- 
ward this end, it is desirable to design the rotor system to 
take maximum advantage of the backup bearings and use 
them as true auxiliary bearings to provide support during 
critical situations in a safe and consistent manner. An im- 
portant concern in this regard is the dynamic behavior of the 
rotor when it comes into contact with the auxiliary bearing. 
If safe and effective operation of the engine is to be ensured 
during these periods, it is essential that designers have a very 
good understanding of the steady-state dynamics of rotor 
systems with clearance effects. 

There are a number of studies in the literature concerned 
with the dynamics of rotors with clearance effects. Ya- 
mamoto (1954) conducted a systematic study of rotor re- 
sponses involving bearing clearance effects. Black (1968) 
studied the rotor/stator interaction with a clearance. He 
concluded that rotor/stator interactions may occur in a vari- 
ety of forms and circumstances, including jump phenomena. 
Ehrich (1966) reported the first identification of a second 
order subharmonic vibration phenomenon in a rotor system 
associated with bearing clearance (1966). Bently (1974) pub- 
lished experimental observations of second and third order 
subharmonic vibration in a rotor system. Later, Muszynska 
(1984) cited the occurrence of second, third, and fourth order 
subharmonic responses in a rotor rubbing case and Ehrich 
(1988 and 1991) observed eighth and ninth order subhar- 
monic vibration as well as chaotic vibration in a high speed 
turbomachine. Childs (1979 and 1982) published two papers 
to explain the mechanism for the second and third order 
subharmonic responses noted above. He stated, with great 
insight, that motion due to nonsymmetric clearance effects 
is a fractional-frequency phenomenon.” 

While those studies have greatly enhanced the understand- 
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FIG. 1 DIAGRAM OF THE FEM ROTOR MODEL 

mg of clearance effects on rotor dynamics, a more detailed 
understanding of the dynamical behavior of such systems is 
needed. The perspective of much of this earlier work is that 
the clearance exists as a result of manufacturing error or mis- 
fitting. That is, it is due to an abnormal situation. However, 
in a rotor system fitted with magnetic bearings and auxiliary 
bearings, the clearance becomes a design parameter rather 
than an irregularity. From this point of view, it is impor- 
tant to develop a detailed quantitative understanding of the 
dynamic responses that are to be expected. Such knowledge 
will provide guidelines for the selection of auxiliary bearing 
parameters. 

It seems that there have been little work to date that is 
specifically concerned with auxiliary bearings in magnetic 
bearing supported rotor systems. Two papers that are di- 
rectly related to research on auxiliary bearings were both fo- 
cused on transient responses. Gelin et al. (1990) studied the 
transient dynamic behavior of rotors on auxiliary bearings 
during the coast down. Ishii and Kirk (1991) investigated 
the transient responses of a flexible rotor during the rotor 
drop after the magnetic bearings become inactive. In both 
papers, idealised rotor models are used and it is assumed 
that once the magnetic bearings fail, the torque is cut off 
and consequently the rotor speed approaches zero. 

In this paper, simulation results are presented for a 
complex rotor system supported by auxiliary bearings with 
clearance at each end of the rotor. This work is specifically 
concerned with systems in which the clearances are quite 
small (on the order of a few mils), which is appropriate 
for jet engine applications in which the backup bearing is 
acting to provide rotor support on a consistent basis. The 
influence of rotor imbalance, support stiffness and support 
damping are investigated using direct numerical integration 
of the governing equations of motion and the harmonic 
balance method. Some insights are obtained with regard 
to the frequency and amplitude behavior of the steady-state 
vibration of such a system. 
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coordinate transformation: 



FIG. 2 AUXILIARY BEARING MODEL 
SIMULATION MODEL 

The rotor is modelled using free-free normal mode shapes 
and natural frequencies obtained through finite element 
analysis. The model data is representative for the rotor 
of a jet engine. Fig. 1 shows a schematic diagram of the 
FEM rotor model. Parametric information about the model 
is listed in Table 1. The torsional motion of the shaft is not 
considered in this paper. Using state space representation 
and modal coordinates, the equations of motion for the rotor 
are expressed as 


N 

XRk = Qxi, 
1 = 1 
N 

<&,t Q ri , 

1 = 1 


(k= NB1,NB2) 


The equations of motion for the auxiliary bearings are 
derived using the model shown in Fig. 2 


M BkXsk 4 Csh^Bk 4* X BkX Bk 

= F nk cos a* - F tk sin a k + M B k g, (2 .a) 
^BkYsk 4 ^Bk^Bk 4 I^Bk^^Bk 

= F n k sin c*k 4 F t k cos a kt (2.6) 

J Bk'&Bk 4 CBrpk^pBk 

Ftk^Bk F n k Rmk t (2.c) 


where 


a* = tan 1 


Yflk — Vgjfc 

XRk — Xsk 


(k = NB1.NB2) 


At this point, the rotor and the back-up bearings appear 
to be uncoupled. However, the force vectors F* and Fy 
on the right hand sides of equations (1) are partially due to 
rotor/ auxiliary bearing interaction. In fact, we have 


where 


Qx42£w a Qx +QI7Qy 4u> 2 Qx 
4 2Q<w n Qy = * T F X , 
QK+2Cu/ n Qy - QTQ X 

-2 QCWnQx =^ r Fy, 


Fx = {Fxi,Fx2,...,Fx m } 1 , 
F Y = 

Qx =^- 1 X /l , 

Qy =$■% 


(la) 

Fxk = 

(1-6) 

Fyk = 


- F n k cos ock 4- F tk sin a k 
4* Mk g 4 M k e Q 2 cos (fit), 

— F n k sina k — F t jt cosa k 
4 Mk e ft 2 sin (ftf). 


The rotor /bearing interaction is represented with the normal 
force F n k 


where 


Fnk = 


( Kc&k, 

10. 


A* <0, 

At >0, 


(3 -a) 


with 

Xft = {-Yftl, Xrj, . . . , ATflm} -1 , 

Y « = 

(m = total number of nodes) 

The physical displacements of the rotor at the two aux- 
iliary bearing locations can be obtained using the following 


A* — (Xft k — Xbi fc) cos ot k 4 (Yfik — Y Bk ) sin a k — S k 

and the Coulomb friction force F tk . As long as there exists 
slip at the contact point, the friction force obeys 

F ik = fiFnk . (3.6) 

However, when there is no slip at the contact point, the 
friction forces are solved from equations (1) and (2) using 
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the kinematic constraint that the circumferential velocities 
of the rotor and the inner-race of the back-up bearing at 
the contact point equal to each other. At the same time, if 
this solved friction force exceeds the maximum static friction 
force (= /ij F nig ), equation (3.b) applys again. 

DISCUSSION OF RESULTS 

The rotor is modeled with 34 stations (as shown in Fig. 1) 
and the first four modes (two rigid body and two flexible 
modes) are included in the simulation model. The two 
auxiliary bearings are located at nodes 3 and 33, respectively. 
This arrangement is taken to represent one of the most 
technically feasible configurations in that it greatly simplifies 
bearing maintenance. It is assumed that the two auxiliary 
bearings are identical in terms of stiffness, damping and 
friction characteristics. Some nominal system parameters 
used for the simulation study are Ac=2.855e+6, £=0.03, 
Rmk=l-1 RBk} /*j= 0.5 and /i^=0.002. To avoid excessive 
cluttering of plots, all the results that are presented in this 
paper correspond to node 3, the location of bearing #1. 

Since the total system which includes two bearings and 
associated friction forces as well as the inner-race motions 
is rather complicated and requires considerable amount of 
computer time for the solutions to converge, the friction 
effect is examined first to see if the model can be further 
simplified. It turns out that the steady-state results obtained 
with and without friction are virtually identical. Even the 
differences in transient responses are quite small, as can be 
seen in Fig. 3(a) and 3(b). The only remarkable effect is 
on the transient responses of the inner-race as shown in 
Fig. 3(c). This observation is confirmed by numerous runs 
using different system parameters and rotor speeds. The 
lack of significance of friction may be attributed to several 
factors. First, the inertia of the inner-race is quite small in 
comparison to the rotor mass, so vibration of the bearing 
has little influence on the rotor vibration. Second, ball 
bearings exhibit quite negligible torsional resistance under 
normal conditions. As a result, the terms that are related 
to the friction forces and the rotational motion of the inner 
race are not included in the simulations that are discussed in 
the following paragraphs. 

The steady-state response characteristics of the system 
are obtained through numerical integration of the simplified 
version of governing equations (1) and (2). Near-zero initial 
conditions are used, simulating situations where the AMBs 
are functioning properly prior to a system failure. Multiple 
solutions with other initial conditions are not sought at 
present. 

It is well known from linear analyses that imbalance 
greatly affects the steady state vibration amplitudes of a 
rotor system. However, it is observed from the current 
work that imbalance may also influence frequency content of 


the rotor responses quite dramatically at certain operating 
speeds. A typical case with such an imbalance effect is shown 
in Figs. 4 and 5, where orbits and corresponding frequency 
spectra of the rotor for different values of imbalance at the 
speed of 0 = 1000 are plotted. For this particular case, 
there exist eight ranges of imbalance values that result in 
eight different types of rotor responses. 

For e < 0.0009, the rotor rotates near the bottom of 
auxiliary bearings and the responses are predominantly syn- 
chronous (Figs. 4(a)-4(b) and 5(a)-5(b)). As imbalance in- 
creases, the 2f2 superharmonic component approaches the 
magnitude order of the synchronous component (Figs. 4(b) 
and 5(b)). However, the responses are of small amplitude. 
For 0.0010 < € < 0.0013, the responses are dominated by 
0/2 subharmonic components (Figs. 4(c) and 5(c)). In other 
words, the amplitude of 0/2 component is greater than that 
of the synchronous. For 0.0014 < e < 0.0016, the 0/2 
subharmonics disappear and the 0/3 subharmonics become 
dominant (Figs. 4(d) and 5(d)). So far, the overall amplitude 
of the responses are not large, the rotor just bounces near the 
bottom of the auxiliary bearings. For 0.0017 < e < 0.0027, 
the orbits become chaotic- looking (Fig. 4(e)) and the spec- 
trum contains a lot of noise (Fig. 5(e)). In this range of 
imbalance, the rotor changes from bouncing near the bot- 
tom to bouncing around the full clearance of the bearing as 
imbalance increases. In the middle of this transition range, 
true chaos is observed. The Poincare map shown in Fig. 6(a) 
and the frequency spectrum shown in Fig. 6(b) demonstrate 
that the response has all the characteristics of a chaotic 
phenomenon. It should be noted that even though the or- 
bits are chaotic looking, the amplitudes are not the largest 
among all the cases for this particular parametric configu- 
ration. For 0.0028 < e < 0.0034, the orbits orbits are no 
longer chaotic-looking (Fig. 4(f)). The spectrum shows they 
are Q/5 subharmonic responses (Fig. 5(f)). Notice the am- 
plitudes are the largest for this parametric configuration. For 
e = 0.0035, the amplitude suddenly becomes smaller even 
though the imbalance has become larger (Fig. 4(g)). And 
the rotor bounces near the bottom of the auxiliary bearings 
again. The frequency spectrum shows it is 0/8 subharmonic 
response (Fig. 5(g)). For 0.0036 < e < 0.0042, the orbits 
become chaotic-looking again (Fig. 4(h)). But the frequency 
spectra are very similar to the subharmonic cases (Fig. 5(h)), 
only with some discrete noise. Finally, for e > 0.0043, the re- 
sponses become predominantly synchronous again (Figs. 4(i)- 
4(j) and 5(i)-5(j)). But this time as imbalance increases, the 
20 superharmonic component become smaller and smaller 
(Figs. 4(j) and 5(j)). 

Examining all the orbits in terms of amplitudes as imbal- 
ance increases, we can see the char act eristics of a jump-type 
phenomenon (Cunningham, 1958). The jump-down takes 
place around 0.0034 < e < 0035 where the rotor jumps 
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from full-clearance bouncing to near— bottom bouncing. Fur- 
ther investigation is needed to better understand this type 
of change. 

Imbalance responses at some other operating speeds and 
for other parametric configurations exhibit similar changes 
as imbalance varies, though the imbalance ranges and cor- 
responding response types may not be as well defined as 
in the above cases. In fact, subharmonic responses from 
f2/2 through fi/ 10 are all observed except 17/9 as shown 
in Figs. 7. Surprisingly, those subharmonics are not di- 
rectly related to the system’s natural frequencies as were the 
cases with other researchers’ findings (such as Ehrich, 1988). 
Moreover, several types of subharmonic responses may occur 
at a single operating speed. It should be noted that Chen et 
al. (1993) also reported occurrence of three stable subhar- 
monic responses at a single rotor speed in a SFD supported 
rotor system but did not provide any explanation for their 
findings. In their case, even the imbalance did not vary. Ap- 
parently, further research is needed to find a mechanism to 
explain this multi-subharmonic vibration phenomenon. On 
the other hand, it should be pointed out that these subhar- 
monic responses are not typical cases. While some of them 
are observed to exist within a certain range of parameters, 
the majority of them occur only for some specific parametric 
configurations. 

Due to space limitations, results for other parametric 
configurations are not systematically plotted. A general 
summation of the observations is presented instead. A 
common feature among all the responses is that for very 
small imbalance, the responses are always synchronous. 
The imbalance range that result in synchronous dominated 
responses depend on several system parameters. For small 
back-up bearing stiffness (such as Kb = 0.213e+6) and 
normal damping ( Cb = 157.0), the responses are almost 
always synchronous. Only fi/2 subharmonic are observed 
at a few operating speeds with a very narrow range of 
imbalance. It should be noted that even though a lower 
Kb may leads to a better system response, it may also 
fail to protect the magnetic bearings due to the fact that 
it could result in a larger rotor orbit-center offset. However, 
the dramatic response changes discussed above may occur 
again if the damping becomes small (such as Cb = 57.0) 
even though the stiffness still remains small. On the other 
hand, increasing the damping Cb alone may not be able 
to eliminate those dramatic changes. It is observed that 
those changes can still occur for Cb being as large as 
700.0. Reducing the size of clearance <5 may not eliminate 
the response changes at certain speeds. But it can narrow 
the operating speed range where those changes occur. For 
example, response changes are eliminated for 0 > 1500 
when 6 is reduced from 0.002 to 0.001 with all other 
parameters remaining the same, but response changes still 


occur for fl < 1400. 

It is obvious that rotor responses involving nonsymmelric 
bearing clearance effect are very complex problems and 
numerical integration alone is not a sufficient tool to obtain 
a global picture of the system responses. The harmonic 
balance method is then used for the investigation of global 
system behavior. However, it is only attempted for situations 
with very small imbalance values. In addition, only the lQ 
harmonic is considered. The complex frequency contents 
associated with medium and large imbalance values makes 
it a formidable task to apply the harmonic balance method 
for other cases. Nevertheless, some useful information can be 
drawn from these results. After ail, an adequately balanced 
rotor system should have very small imbalance under normal 
conditions. 

Fig. 8(a) shows that nonsymmetric clearance effect is 
equivalent to asymmetric stiffness effect with regards to 
critical speeds. The clearance actually splits the first critical 
speed into two pseudo-critical speeds. In the X direction, 
the gravity force tends to keep the rotor in contact with 
the bearing at low operating speed. Thus, the apparent 
stiffness is almost the same as Kb and the pseudo-critical 
speed is nearly the same as the critical speed for the linear 
case (£ — 0). But in the Y direction, the clearance 
results in a lower apparent stiffness and, consequently, an 
additional lower-value pseudo-critical speed. It is seen that 
several higher order additional pseudo-critical speeds are 
created in the operating speed range in addition to the 
1st additional pseudo-critical speed. It is noted that the 
response in the X direction also departs from the linear 
case at high operating speed. This is because the imbalance 
force becomes dominant at high rotor speed which in torn 
makes the gravity force less significant and the clearance 
effect more important. Fig. 8(b) shows that changing the 
auxiliary bearing stiffness has little effect on the pseudo- 
critical speeds of the system. However, for a larger value 
of imbalance, a higher Kb does leads to a greater tendency 
of double-valued responses. In each direction, for either a 
stiffness increase or an imbalance increase, the 1st pseudo- 
critical peak tend to become a hard-spring type jump aud 
the 2nd one tend to develop into a soft-spring type jump, 
with the tendency decreasing as the pseudo— critical’s order 
increases. 

Fig. 9(a) shows double-valued responses in the Y direction 
for four different values of clearance. It is seen that a larger 
clearance results in wider rotor speed range of double-valued 
responses. It is also observed that as clearance incresses, 
the apparent stiffness decreases and the first pseudo-critical 
speed shifts to a lower value. Fig. 9(b) shows the double- 
valued responses in the X direction. Even though the jumps 
themselves are smaller in magnitude, they are more obvious 
in trend. Notice how little the change is for the first pseudo- 
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critical speed in the X direction. 

Fig. 10(a) shows the influence of auxiliary bearing damping 
on the double- valued responses in the Y direction. It is 
observed that the damping has to be quite large to eliminate 
the double-valued responses associated with the first pseudo- 
critical speed. Fig. 10(b) shows the influence of C B on the 
double-valued responses in the X direction. In both figures, 
it should be noted that as C B decreases, the second pseudo- 
critical speed peak will develop into a soft-spring type jump 
and the third pseudo-critical speed peak will evolve into a 
hard-spring type jump. 

The system behaviors for higher operating speed range are 
not shown in Figs. (9) and (10) so that the jump phenomena 
can be more clearly illustrated. It is also because that the 
system’s responses at high operating speed range with the 
same parameters are more or less regular, in other words, 
mainly amplitude changes. 

CONCLUSIONS 

As a summary of the results discussed above, the following 
conclusions can be drawn: 

1* Ih* balance may change the rotor responses dramatically in 
terms of frequency contents at certain operating speeds, 
especially under conditions of large clearance, high bearing 
stiffness and low bearing damping. With imbalance 
changing, as many as eight different types of responses 
may occur for a particular parametric configuration at a 
single operating speed. 

2. Subharmonic responses of second order through tenth 
order are all observed except for the ninth order case. 
However, the majority of them are not typical cases, 
and were observed only for quite particular parametric 
configurations. 

3. Chaotic phenomenon is observed to occur occasionally. 
However, the amplitudes associated with such motion are 
not among the largest. 

4. Nonsymmetric clearance effects influence the critical speeds 
in a manner similar to asymmetric stiffness effects. 

5. Double-valued responses in the form of both hard-spring 
type jump and soft-spring type jump are observed to be 
possible at low operating speeds with low auxiliary bearing 
damping or high imbalance. With large clearances or high 
bearing stiffness, the jump phenomenon may occur for 
even relatively small imbalances. 

6 . The effect of friction between the shaft and the inner race 
of a roiling element auxiliary bearing on the dynamics of 
the rotor is quite small and can reasonably be neglected 
for steady-state analyses. 


discussions and practical advice. 
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TABLE 1. FEM MODEL DATA 


Station 

(in.) 

O.D. 

(in.) 

ID. 

(in.) 

WT. 

(lb) 

Polar Mom. 
(Ib.in.s 2 ) 

5.400 

2.800 

2.000 

0.000 

0.000 

5.450 

2.800 

2.000 

5.000 

5.000 

6.500 

2.800 

2.000 

0.000 

0.000 

6.505 

7.300 

6.950 

0.100 

0.100 

7.100 

7.300 

6.950 

0.000 

0.000 

8.670 

7.300 

7.120 

30.400 

425.000 

10.900 

8.400 

8.220 

0.000 

0.000 

13.120 

9.300 

9.120 

0.000 

0.000 

15.310 

10.100 

9.920 

20.300 

538.000 

17.280 

10.700 

10.550 

9.150 

302.000 

18.580 

11.200 

11.050 

0.000 

0.000 

19.880 

11.550 

11.400 

0.000 

0.000 

21.120 

11.950 

11.800 

0.000 

0.000 

22.380 

12.150 

12.000 

41.900 

1645.000 

23.670 

12.400 

12.250 

0.000 

0.000 

24.900 

12.500 

12.350 

0.000 

0.000 

26.120 

12.500 

12.350 

0.000 

0.000 

27.400 

12.500 

12.350 

22.550 

539.000 

27.750 

12.500 

2.500 

1.027 

0.000 

27.800 

4.500 

2.500 

7.733 

0.000 

32.000 

4.500 

2.500 

25.022 

49.761 

37.000 

4.500 

2.500 

17.961 

0.000 

43.000 

4.500 

2.500 

17.961 

0.000 

48.000 

4.500 

2.500 

16.328 

0.000 

53.000 

4.500 

2.500 

16.328 

0.000 

58.000 

4.500 

2.500 

14.614 

0.000 

61.950 

4.500 

2.500 

6.531 

0.000 

62.000 

6.000 

2.500 

0.954 

0.000 

62.405 

6.000 

5.670 

34.935 

1216.500 

64.600 

6.000 

5.670 

0.000 

0.000 

67.000 

6.000 

5.670 

117.000 

3750.500 

69.500 

6.000 

5.670 

0.000 

0.000 

69.505 

2.850 

1.250 

0.000 

0.000 

72.500 

2.850 

1.250 

10.197 

12.000 

73.600 

2.850 

1.250 

0.000 

0.000 


Modulus of elasticity: E=2.80e+7 psi 
Shear Modulus of elasticity: G=1.08e+7 psi 
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FIG. 3 EFFECTS OF BEARING FRICTION (C B =157.0, 
/C B =0.313e+6, 0=1400, <5=0.002, e=0.0002) 
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FIG. 7 SUBHARMONIC RESPONSES (ORBITS) 
(/C B =0.313e+6) 
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FIG. 8 INFLUENCE ON CRITICAL SPEEDS 
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(e=0.0001, A' s =0.313e+6, C s =300.0) 
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FIG. 10 EFFECTS OF BACK-UP BEARING DAMPING 
(6=0.002, e=0.0001, A B =0.313e+6) 
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FIG. 6 CHAOTIC RESPONSES (e = 0.0023, 
6=0.002, Afi=0.313e+6, C B = 157.0) 











